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HIGHLIGHTS 


•  DCFC  developed  here  uses  a  porous  alumina  tube  to  contain  carbon/carbonate  slurry. 

•  Bubbling  in  carbon/carbonate  slurry  impact  is  dictated  by  the  C/02-  ratio. 

•  Bubbling  inhibits  partial  electrochemical  oxidation  when  the  C/02-  ratio  is  low. 

•  Bubbling  increases  total  efficiency  from  52  %  to  64  %  for  WC/ca rbonate  =  1.0  wt%. 
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This  study  aims  to  enhance  the  direct  carbon  fuel  cell  (DCFC)  efficiency  by  bubbling  Ar  gas  in  a  carbon/ 
carbonate  slurry.  The  current  discharge  curve  and  CO  and  CO2  production  rates  were  measured  during 
discharge  at  20  mA  cm-2.  When  carbon  content  in  carbonates  (Wc/carbonate)  was  1.0  wt%,  Ar  bubbling 
decreased  the  CO/CO2  production  ratio  from  0.16  to  0.051,  and  increased  the  CO2  production  rate,  which 
approached  the  theoretical  value  for  complete  electrochemical  oxidation  of  carbon.  Moreover,  bubbling 
increased  the  cell  voltage  by  about  0.1  V  at  steady  state.  At  Wc/carbonate  =  3.0  wt%,  the  CO/CO2  production 
ratio  without  bubbling  increased  significantly  to  8.5,  indicating  that  partial  electrochemical  oxidation 
became  dominant.  The  C/O2-  ratio  is  key  for  explaining  the  electrochemical  oxidation  of  the  carbon. 
When  the  C/O2-  ratio  was  low,  such  as  at  Wc/carbonate  =  1.0  wt%,  complete  electrochemical  oxidation 
became  dominant,  and  bubbling  enhanced  the  complete  electrochemical  oxidation.  When  the  C/O2- 
ratio  was  high,  such  as  at  Wc/carbonate  =  3.0  wt%,  partial  electrochemical  oxidation  became  dominant,  and 
bubbling  did  not  inhibit  it.  For  Wc/carbonate  =  1-0  wt%,  bubbling  increased  both  coulombic  and  voltage 
efficiencies,  resulting  in  an  increase  in  total  efficiency  from  52%  to  64%. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  able  to  efficiently  convert  the  chemical  energy 
from  fuels  into  electrical  energy.  However,  conventional  fuel  cells 
like  proton-exchange  fuel  cells  (PEFCs)  and  solid  oxide  fuel  cells 
(SOFCs)  use  gaseous  fuels.  Conversion  of  energy  from  solid  carbon 
fuels  such  as  coal  or  biomass  chars  in  these  cells  requires  gasifi¬ 
cation.  In  contrast,  direct  carbon  fuel  cells  (DCFCs)  convert  the 
chemical  energy  in  carbon  directly  into  electricity,  without  the 
need  for  gasification.  The  solid  carbon  fuels  used  in  DCFCs  can  be 
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produced  from  sources  such  as  coal,  natural  gas,  and  biomass.  Fine 
carbon  particles  are  electrochemically  oxidized  at  high  tempera¬ 
tures  (C  +  O2  -►  CO2). 

Although  variety  of  DCFC  designs  have  been  tested,  and  some 
articles  related  to  DCFCs  have  recently  been  published  [1—3  ,  this 
technology  is  still  in  the  preliminary  stage  of  development  [1  . 
There  have  been  various  attempts  to  convert  carbon  materials 
directly  into  electricity  in  DCFCs  using  electrolytes,  such  as  molten 
carbonates  [4-7  ,  molten  hydroxides  [8  ,  and  YSZ-based  solid 
electrolytes  [9,10  .  The  use  of  molten  carbonate  has  a  number  of 
advantages:  it  has  good  long-term  stability  in  CO2,  can  catalyze 
carbon  oxidation,  and  has  a  high  ionic  conductivity  1  .  When 
carbon  particles  such  as  pulverized  chars  are  used  in  DCFCs,  carbon/ 
carbonate  slurry  is  an  appropriate  fuel  for  improving  contact  with 
metal  current  collector  or  anode. 
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The  anode  (Rl),  cathode  (R2),  and  overall  reaction  (R3)  in  this 
system  are  expressed  as  follows: 

C  +  2CO§“— 3C02  +4e~  (Rl) 

02  +  2C02+4e-^2C0^-  (R2) 

C  +  02  C02  (R3) 

The  desired  reaction  is  complete  oxidation,  which  releases  four 
electrons  per  carbon,  as  given  by  Rl.  However,  several  studies  have 
shown  a  possible  formation  of  CO  at  the  anode  compartment  from 
R4,  R5,  or  even  from  the  Boudouard  reaction  (R6)  [4,6,11-14  . 

C  +  C03~^C0  +  C02  +  2e~  (R4) 

C+  l/2CO^“->3/2CO  +  e-  (R5) 

C  +  C02  ->•  2 CO  (R6) 

Two  electrons  are  released  per  carbon  atom  in  R4,  whereas  only 
one  electron  is  released  per  carbon  atom  in  R5.  However,  this  aspect 
remains  speculative  because  gas  analysis  was  not  always  provided 
or  performed  in  the  relevant  research  studies  [11  . 

The  total  electrochemical  efficiency  of  a  DCFC  is  determined  by  its 
voltage  efficiency  (the  ratio  of  cell  voltage  to  theoretical  voltage)  and 
coulombic  efficiency  (electron  yield  (n)  per  carbon  atom,  divided  by 
4)  [15  .  Cooper  and  Selman  derived  a  coulombic  efficiency  formula  in 
DCFC  from  two  reactions  (Rl  and  R4)  16],  and  a  recent  study  has 
employed  coulombic  efficiency  data  to  evaluate  DCFC  performance 
[17];  however,  little  effort  has  been  made  to  enhance  coulombic  ef¬ 
ficiency  by  inhibiting  partial  electrochemical  oxidations  (R4  and  R5). 

This  study  aims  to  enhance  DCFC  efficiency  by  bubbling  an  inert 
gas  in  a  carbon/carbonate  slurry.  Inert  gas  bubbling  can  be  easily 
applied  to  stir  highly  corrosive  molten  carbonates.  The  current 
discharge  curve  and  CO  and  CO2  production  rates  were  measured, 
and  partial  electrochemical  oxidations  were  studied. 

2.  Experimental  section 

2.1.  Carbon  characterizations 

Commercial  activated  carbon  (activated  charcoal,  Kantokagaku) 
was  used  in  this  study.  Elemental  composition  of  activated  carbon 
was  determined  using  CHN  elemental  analyzer  (Micro  corder  JM10, 
J-Science  Lab  Co.  Ltd.)  and  sulfur  analyzer  (HSU-20,  J-Science  Lab 
Co.  Ltd.).  Characterization  of  the  activated  carbon  was  performed  by 
X-ray  diffraction  (X'Pert  MPD,  Philips)  with  Cu  Ka  radiation  at  40  kV 
and  30  mA.  The  size  distribution  of  activated  carbon  was  measured 
by  laser  diffraction  (LMS-2000e,  Seishin  Enterprise  Co.,  Ltd.). 

2.2.  DCFC  experiments 

Lig.  1  shows  a  schematic  diagram  of  a  DCLC  with  a  molten  car¬ 
bonate  electrolyte,  which  was  developed  in  our  laboratory.  In  this 
study,  the  working  electrode  (WE),  counter  electrode  (CE),  and 
reference  electrode  (RE)  were  made  from  gold  sheet.  The  area  in 
contact  with  the  carbon  was  1.0  cm2.  The  RE  and  CE  were  in  alumina 
tubes,  and  the  WE  was  in  a  porous  alumina  tube  (for  separation  of 
carbon  particles,  as  described  later).  Each  gold  sheet  was  spot  welded 
to  the  gold  wire  and  extended  to  the  other  end  of  the  tube  to  provide  a 
connection  to  the  potentiostat/galvanostat  (HAL-3001,  Hokuto 
Denko)  for  cell  parameter  measurements.  On  the  basis  of  cell  mea¬ 
surements,  current  density  was  set  at  20  mA  cm'2  during  discharge. 


Carbonate  easily  corrodes  quartz;  therefore,  a  high  purity  alumina 
crucible  was  used  as  a  holder  for  carbonate.  In  order  to  keep  the  entire 
cell  in  a  gastight  environment,  the  alumina  crucible  was  placed  at  the 
bottom  of  the  quartz  reactor  with  O-rings  and  silicone  plugs.  The  gas 
feed  tubes  for  the  counter  electrode  (CE)  and  the  reference  electrode 
(RE),  and  the  porous  tube  for  the  working  electrode  (WE)  were  made 
from  high  purity  alumina,  and  these  were  supported  by  silicone  plugs 
attached  to  the  top  of  the  quartz  reactor  as  a  gas  seal.  The  alumina 
crucible  contained  243  g  of  dry  ternary  carbonate  powder 
(Li2C03:Na2C03:I<2C03  =  1:4:4  (weight  ratio))  without  carbon,  and 
the  mixed  carbon/carbonate  powder  was  contained  in  the  porous 
alumina  tube  having  an  average  pore  diameter  of  95.6  nm.  Therefore, 
carbon  particles  larger  than  this  size  did  not  pass  through  the  porous 
alumina  tube.  This  allowed  contact  between  the  molten  carbonate 
electrolyte  and  WE,  CE,  and  RE.  In  the  porous  alumina  tube,  the  mixed 
carbon/carbonate  powder  consisted  of  11.3  g  of  dry  ternary  carbonate 
powder  (Li2C03:Na2C03:I<2C03  =  1 :4:4  (weight  ratio))  and  0.113  g  or 
0.339  g  of  activated  carbon  (i.e.,  1.0  wt%  or  3.0  wt%  of  the  carbon  fuel  in 
the  carbonates,  defined  as  Wc/Carbonate)-  Although  the  single  carbonate 
has  relatively  high  melting  point:  L^CCU  (999  K),  Na2C03  (1131  K),  and 
K2CO3  (1172  K),  the  eutectic  carbonate  has  much  lower  melting  point; 
therefore,  ternary  eutectic  carbonate  has  been  generally  used  in  the 
DCFC  with  a  molten  carbonate  electrolyte  [7]. 

After  the  fuel  cell  was  assembled  and  sealed  to  be  gastight,  the 
reactor  was  heated  using  an  electric  furnace  to  1073  K.  Ar 
(100  ml  min-1)  was  introduced  into  the  WE  compartment  during 
heating.  After  the  carbonate  powder  melted,  alumina  tubes  with  CE 
and  RE  were  inserted  into  the  molten  carbonates  in  the  alumina 
crucible.  Then,  02  (50  ml  min-1)  and  C02  (100  ml  min-1)  gases 
were  introduced.  The  DCFC  developed  here  was  similar  to  that 
developed  in  previous  studies  [4,6,12,14] ;  but,  there  were  some 
differences.  In  the  DCFCs  used  in  previous  studies,  carbon/carbon¬ 
ate  slurry  contacted  the  WE  that  was  contained  in  an  alumina 
crucible,  and  pure  carbonates  electrolyte  contacted  the  CE  and  RE 
that  was  contained  in  an  alumina  tube  with  holes.  In  contrast,  here 
the  carbon/carbonate  slurry  was  contained  in  the  porous  alumina 
tube,  and  pure  carbonates  were  contained  in  the  alumina  crucible. 
The  positions  of  the  carbon/carbonate  slurry  and  the  carbonate  in 
the  DCFCs  in  this  study  were  opposite  to  the  positions  in  DCFCs  in 
previous  studies  4,6,12,14  .  Containing  the  carbon/carbonate  slurry 
in  the  porous  alumina  tube  is  advantageous  because  the  slurry  is 
easily  exchanged  by  using  new  porous  alumina  tube  with  the  slurry 
when  ash  accumulates  over  time  with  the  use  of  coal  chars.  In  this 
study,  activated  carbon  was  used  as  fuel.  However,  coal  chars  were 
used  as  well.  Build-up  of  ash  in  the  molten  carbonate  electrolyte  is 
one  of  the  main  barriers  in  DCFC  development  [5].  The  DCFC 
developed  in  this  study  is  of  a  cartridge  type,  and  the  cell  stack  is 
packaged  by  using  several  anode  compartments  consisting  of 
porous  alumina  tubes  and  carbon/carbonate  slurry. 

This  DCFC  comprised  WE,  CE,  RE,  and  reactor  compartments. 
Herein,  a  reactor  compartment  is  a  space  in  the  quartz  reactor, 
excluding  the  WE,  CE,  and  RE  compartments.  Anode  off  gases  from  WE 
were  able  to  pass  through  the  porous  alumina  tube,  and  diffuse  to  the 
reactor  compartments.  Therefore,  after  reaching  steady  state,  CO2  and 
CO  concentrations  in  anode-off  gases  from  WE  and  reactor  compart¬ 
ments  (XC02,WE,  XC02, reactor,  Xco,we,  XCo, reactor)  were  measured  with  a 
flame  ionization  detector  (FID)  equipped  with  a  Porapack-Q  column 
and  methanizer.  Outlet  gas  flow  rates  (Qwe,  and  Qreactor)  were  also 
measured,  using  a  wet  gas  meter  (W-NK-1,  Shinagawa  Co.,  Ltd). 

CO2  and  CO  flow  rates  (Qco2  and  Qco)  are  given  as  follows: 

Qco2=E(Xco2,iQi  Inc  discharge)  (i  =  WE>  reactor  compartment) 
i 


(i) 
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Fig.  1.  Schematic  diagram  of  DCFC  with  a  molten  carbonate  electrolyte  using  porous  alumina  tube  in  anode  compartment. 
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alumina 

tube 


Qco  =  Y,  (xco,iQj|no  discharge)  ('  =  WE>  reactor  compartment) 
i 

(2) 

C02  and  CO  production  rates  (Qc02, product  and  Qco, product)  are 
determined  by  the  sum  of  the  differences  between  C02  and  CO  flow 
rates  during  discharge  (at  20  mA  cnrT2)  and  those  at  the  OCV  (open 
circuit  voltage)  condition,  given  as: 


0cO2,  product  —  (Xc°2  4  Qj I  discharge  ^C02  4  Q* I  no  discharge) 

i 

(i  =  WE,  reactor  compartment) 

Qco, product  —  (j^CO, iQi  (discharge  ^CO.iQi  |no  discharge) 

i 

(i  =  WE,  reactor  compartment) 

The  C02  production  rate  at  the  time  discharge  is  processed 
through  R1  and  at  constant  is  calculated  as  follows: 


(3) 

(4] 


0cO2, product  —  (3) 

where  F  is  the  Faraday  constant  (96,450  C  mol”1),  Vm  is  the  molar 
volume  of  C02  at  standard  temperature  and  pressure,  and  I  is  the 
constant  current  flow.  Unit  of  Ampere  corresponds  to  coulomb  per 
second  (C  s”1).  The  C02  production  rate  is  0.209  ml  min”1  when  I  is 
20  mA.  Anode-off  gases  were  measured  before  discharge,  con¬ 
firming  that  02  concentration  was  negligible  in  the  anode  and 
reactor  compartments. 

2.3.  Ar  bubbling  in  anode  compartment 

Fig.  2  shows  a  schematic  diagram  of  the  anode  compartment 
with  (a)  and  without  (b)  bubbling.  Ar  gas  for  bubbling  was  released 
from  two  holes  in  a  gas  line  tube.  The  diameters  of  the  holes  were 
2  mm.  When  the  carbon/carbonate  slurry  was  used  as  fuel  in  DCFCs, 
fuel  stirring  played  an  important  role  in  determining  the  electro¬ 
chemical  performance  [4,6].  In  this  study,  fuel  was  stirred  using  Ar 
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(b)  Without  bubbling 


Fig.  2.  Schematic  diagram  of  anode  compartment  with  or  without  bubbling. 


gas  bubbling  with  a  flow  rate  of  50  ml  min^1.  For  stirring  molten 
carbonates,  high  purity  aluminum  or  other  materials  resistant  to 
highly  corrosive  molten  carbonates  are  required.  In  contrast,  inert 
gas  bubbling  is  easier  to  use  with  molten  carbonates.  Bubbling  is 
supposed  to  stir  the  carbon/carbonate  slurry,  and  provide  carbons 
and  ions  to  the  anode  surface.  Unlike  mechanical  stirring,  bubbling 
is  also  expected  to  immediately  remove  gases  (CO2  and  CO) 
resulting  from  the  electrochemical  reaction. 

3.  Results  and  discussion 

3.2.  Carbon  characterizations 

Table  1  shows  an  elemental  analysis  (dry-base)  and  the  ash 
content  of  samples. 

Fig.  3  shows  the  XRD  pattern  of  the  activated  carbon.  Graphite 
carbon  was  also  analyzed  for  reference,  and  an  intense  peak  corre¬ 
sponding  to  the  (002)  graphitic  basal  plane  reflection  was  observed  at 
20  =  26.5.  These  peaks  are  relatively  low,  indicating  disordered 
structure  than  graphite  carbon,  which  is  consistent  with  previous 


studies  [6,14].  The  average  size  of  carbon  crystalline  was  calculated 
from  the  Debye-Scherer  equation  [6  .  Constants  in  the  Debye-Sch- 
erer  equation  were  0.9  for  Lc  (the  layer  dimension  perpendicular  to 
the  basal  plane)  and  1.84  for  La  (the  layer  dimension  parallel  to  the 
basal  plane).  The  crystalline  parameters,  d0 02  (interplanar  distance), 
Lc,  and  La  are  shown  in  Table  2.  The  relationship  between  carbon 
properties  and  OCV  in  the  DCFC  will  be  discussed  later. 

Fig.  4  shows  the  size  distribution  of  activated  carbon.  The  average 
diameter  of  activated  carbon  was  22.4  pm.  Thus,  the  activated  carbon 
did  not  pass  through  the  porous  alumina  tube  used  in  the  DCFC  cell. 

3.2.  Anode  offgas  characteristics 

Fig.  5  shows  measured  CO2  and  CO  flow  rates  without  bubbling 
at  the  OCV  condition  (Jd  =  0  mA  cnrT2)  and  Jd  =  20  mA  cm-2,  and  CO 


Table  1 

Elemental  analysis  of  the  activated  chars  [wt%,  dry]. 


c 

H 

N 

S 

O  (diff) 

Ash 

92.3 

0.17 

3.04 

0.10 

4.39 

0.8 
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Fig.  3.  XRD  patterns  of  activated  and  graphite  carbon. 
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estimation  at  Wc/Carbonate  =  1-0  wt%.  Some  carbonate  was  contained 
in  the  alumina  crucible;  therefore,  CO2  was  formed  at  the  OCV 
condition  as  a  consequence  of  the  carbonate  decompositions  given 

by  R7  [18-20]; 

M2C03  -►  M20  +  C02  (R7) 

where  M  is  K,  Na  or  Li.  CO2  formation  at  the  OCV  condition  also 
indicates  M2O  formation.  Especially,  the  single  Li2C03  starts  to 
decompose  at  lower  temperature  than  the  reaction  temperature  of 
1073  K  20].  L^O  is  an  important  species  to  achieve  a  higher  OCV  in 
DCFCs  as  described  later.  CO  is  also  produced  at  the  OCV  condition. 
In  addition  to  R6,  CO  forms  through  carbothermic  reactions  given  as 
R8  and  R9  [18-20]. 

M2C03  +  C  -►  M20  +  2 CO  (R8) 

M20  +  C  -►  2M  +  CO  (R9) 

CO2  and  CO  flow  rates  increased  at  7d  =  20  mA  cm-2.  The  in¬ 
crease  in  CO2  comes  from  electrochemical  oxidations  given  as  R1 
and  R4.  There  are  two  possible  oxidation  reactions  resulting  in  the 
CO  increase  during  discharge.  One  is  the  Boudouard  reaction,  and 
the  other  is  partial  electrochemical  oxidation.  An  increase  in  CO2 
during  discharge  may  increase  the  rate  of  the  Boudouard  reaction. 
If  this  is  the  case,  then  the  rate  of  increase  of  CO  is  the  same  as  that 
of  CO2  because  the  reaction  rate  of  R6  has  a  linear  relationship  with 
CO2  concentration.  In  addition  to  this,  here,  CO  is  assumed  to  form 
through  R6  at  OCV  condition.  CO  flow  rate  estimated  from  these 
assumptions  is  described  as  “OCV  +  R6”.  However,  the  measured 
CO  flow  rate  at  /  =  20  mA  was  about  70%  higher  than  the  estimated 
value,  suggesting  that  the  dominant  reaction  in  CO  formation 
during  discharge  is  not  the  Boudouard  reaction  but  partial  elec¬ 
trochemical  oxidation. 

Fig.  6  shows  measured  C02  and  CO  flow  rates  without  bubbling 
in  the  OCV  condition  (Jd  =  0  mA  cm-2)  and  Jd  =  20  mA  cnrT2,  and  CO 


Table  2 

Crystalline  parameters  of  activated  chars. 


d002  [nm] 

Lc  [nm] 

La  [nm] 

0.376 

0.987 

3.22 

estimation  at  Wc/carbonate  =  3.0  wt%.  In  the  OCV  condition,  CO 
concentration  at  Wc/C arbonate  =  3.0  wt%  was  almost  three  times  that 
at  Wc/carbonate  =  1.0  wt%.  This  indicates  that  CO  is  produced  through 
R6,  R8,  and  R9,  because  these  reaction  rates  have  linear  relations 
with  carbon  content  in  carbonates.  At  Wc/carbonate  =  3.0  wt%,  the 
discharge  process  does  not  increase  C02  flow  rate,  but  it  signifi¬ 
cantly  increases  CO  flow  rate.  This  is  clear  evidence  that  partial 
electrochemical  oxidation  plays  an  important  role  during  discharge 
at  higher  carbon  content,  as  discussed  in  detail  later  in  this  study. 


3.3.  Effect  of  bubbling  on  current  discharge  curves  and  anode  off 
gases 

Fig.  7  shows  the  effect  of  bubbling  on  constant  current  discharge 
curves  at  WC/c arbonate  =  1-0  wt%.  Voltage  at  t  =  0  min  corresponds  to 
the  OCV  condition.  At  Wq carbonate  =  1.0  wt%,  OCVs  with  and  without 
bubbling  are  1.15  V  and  1.27  V,  respectively  (see  Fig.  7(a)).  The 
standard  potential  of  R1  (C  (graphite)  +  02-^  CO2)  is  1.026  V  as 
evaluated  from  the  Gibbs  energy  change  at  1073  I<  21];  however, 
experimental  OCV  is  generally  higher  than  the  standard  potential  of 
graphite  carbon  as  a  consequence  of  carbon  properties  and  gas 
composition  [5,22,23  .  As  shown  in  Fig.  3,  the  activated  carbon  has 
disordered  structure.  The  highly  pore  crystallized,  lattice  disor¬ 
dered  carbons  are  more  reactive  than  highly  oriented  pyrolytic 
graphite,  resulting  in  higher  OCVs  [22  .  However,  even  if  these 
factors  are  considered;  1.27  V  is  significantly  higher  than  the 
standard  potential.  In  several  studies,  OCV  values  of  1.2-1. 5  V  were 
reported.  Such  values  are  significantly  higher  than  the  thermody¬ 
namically  expected  value  (Eo  =  1.026  V  at  1073  K)  [24-26  .  These 
high  OCV  values  suggest  the  intervening  role  of  CO  and/or  the  alkali 
ions  in  the  molten  carbonate  participating  in  the  anode  reaction 
[25,26].  When  CO  participation  is  considered  as  R10,  the  thermo¬ 
dynamically  expected  value  is  still  lower.  However,  if  Li+  is 
considered  as  the  counter  cation,  net  cell  reactions  are  described  as 
Rll  and  R12,  and  OCVs  are  higher  than  1.026  V: 

CO  +  l/202  ->•  C02  (£°  =  0.98  V  at  1073  K)  (R10) 

C  (s)  +  02  (g)  +  Li20  (cr)  -►  Li2C03  (cr)  (E°  =  1.18  V  at  1073  K)(R11) 
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Fig.  5.  Measured  C02  and  CO  flow  rate  without  anode  bubbling  at  OCV  (/  =  0  mA  cm  2),  /  =  20  mA  cm  2,  and  estimation  (OCV  +  R6)  (lVC/Carbonate  =  1.0  wt%). 


CO  (g)  +  l/202  (g)  +  Li20  (cr)  -  Li2C03  (cr) 

(£°  =  1.30  V  at  1073  I<)  (R12) 

Details  of  Rll  and  R12  are  described  in  Refs.  [25,26].  The  liter¬ 
ature  suggests  that  R12  provides  an  acceptable  explanation  of  high 
OCV  of  1.3  V  [24-26].  As  previously  described,  CO  is  formed  by  R6, 
R8,  and  R9,  and  U2O  is  formed  by  R7  and  R8.  This  also  explains  why 
Ar  bubbling  decreases  OCV.  Ar  bubbling  is  supposed  to  remove 
gases  from  carbon/carbonates  slurry  immediately;  therefore,  Ar 
bubbling  decreases  OCV  to  1.1 -1.2  V  by  removing  CO  before  it 
consumes  through  R12.  Meanwhile,  Ar  bubbling  improves  the 
current  discharge  curve  at  steady  state.  At  the  beginning  of  the 
discharge  process,  the  cell  voltage  rapidly  decreases,  and  then  be¬ 
comes  almost  constant. 

To  calculate  the  theoretical  CO  flow  rate  through  R12,  it  is 
assumed  that  R12  consists  of  following  anode  reaction  (R13). 

CO  +  02~  +  U20  -►  U2C03  +  2e~  (R13) 

Based  on  R13,  the  theoretical  CO  flow  rate  necessary  for 
20  mA  cm-2  discharge  is  0.139  ml  min-1.  However,  CO  flow  rate  is 
0.016  ml  min-1  at  the  OCV  condition,  which  is  much  lower  than  the 
theoretical  value.  Thus,  cell  voltage  decreases  rapidly  when 
bubbling  is  not  used.  Bubbling  increases  the  cell  voltage  at 
approximately  3-5  min  in  current  discharge  curves,  and  cell 
voltage  obtained  with  bubbling  is  higher  than  that  without 
bubbling  by  approximately  0.1  V  at  20  min.  This  is  because  bubbling 
is  supposed  to  stir  the  fuel/carbonate  slurry,  and  as  a  result  it 
provides  both  carbon  and  ions  to  reaction  sites.  In  fact,  the  voltage 
reduction  results  from  diffusion  polarization. 


Fig.  8  shows  the  effect  of  bubbling  on  constant  current  discharge 
curves  at  Wqc arbonate  =  3.0  Wt%.  At  Wqc arbonate  =  3.0  wt%,  the  OCV 
difference  between  with  and  without  bubbling  is  less  than  3%.  It  could 
be  due  to  an  experimental  error.  Finding  here  is  that  bubbling  does  not 
influence  the  current  discharge  curve  in  contrast  to  that  at  Wq 
carbonate  =  1.0  wt%,  and  the  OCV  obtained  is  1.3— 1.4  V.  This  is  because 
CO  concentration  at  Wqc arbonate  =  3.0  wt%  is  about  three  times  higher 
than  that  at  Wqc arbonate  =  1.0  wt%  at  the  OCV  condition.  Although  Ar 
bubbling  promoted  the  removal  of  gases  from  the  carbon/carbonates 
slurry,  complete  removal  of  CO  is  beyond  the  capacity  of  Ar  bubbling. 
Moreover,  CO  concentration  increases  significantly  because  of  partial 
electrochemical  oxidation  during  discharge,  and  it  is  higher  than  the 
0.139  ml  min-1  that  is  necessary  for  a  20  mA  cm-2  discharge. 
Therefore,  R12  plays  an  important  role  in  discharge  at  Wq 
carbonate  =  3.0  wt%,  leading  to  higher  cell  voltage  than  at  Wq 
carbonate  =  1.0  wt%.  Although  it  is  difficult  to  determine  the  minimum 
CO  concentration  that  is  needed  to  have  an  effect  on  the  OCV,  the 
findings  described  herein  provide  acceptable  explanations  for  the 
insignificant  effect  of  Ar  bubbling  on  the  current  discharge  curve  at 
Wc/carbonate  —  3.0  Wt/4. 

Fig.  9  shows  the  effect  of  bubbling  on  measured  CO2  and  CO 
production  rates  during  discharge  at  Wc/Carbonate  =  TO  wt%.  Ar 
bubbling  decreases  the  CO/CO2  production  ratio  from  0.16  to  0.051, 
and  increases  the  CO2  production  rate,  which  approached  the 
theoretical  value  (=0.209  ml  min-1). 

Fig.  10  shows  the  effect  of  bubbling  on  measured  CO2  and  CO 
production  rates  during  discharge  at  Wc/Carbonate  =  3.0  wt%.  The  CO/ 
C02  ratios  with  and  without  bubbling  are  8.5  and  9.0,  respectively. 
In  contrast  to  the  result  at  Wqc arbonate  =  TO  wt%,  bubbling  increases 
the  CO  production  rate  at  Wc/Carbonate  =  3.0  wt%,  because  it  removes 


Fig.  6.  Measured  C02  and  CO  flow  rate  without  anode  bubbling  at  OCV,  20  mA  cm  2  and  estimation  (OCV  +  R6)  (WC/c arbonate  =  3.0  wt%). 


346 


H.  Watanabe  et  al.  /  Journal  of  Power  Sources  273  (2015)  340-350 


O 

> 


1.4 


1.2  - 

o 

1  - 


* 


* 


o 


o; 


^S2o00Sogoo0ooo000o°o0Ooo0oO 

*  **  *  ****  **„ 

7K  ^  7K  ^2  Si/  7R  \i/  ...  \i/  v,.  v..  v/ 


O  °  o  c  ,,  °  °  ° 


*****  ************* 


*  *  *  * 


0.4  - 


0.2  - 


0 


without  bubbling 
with  bubbling 


* 


o 


o 


0 


10 


15 


20 


Elapsed  time  [min] 


C  +  CO2-  -^2C0  +  02~  (R15) 

Electrochemical  reaction  mechanisms  after  02~  formation  are 
summarized  in  previous  studies  [5,29].  Briefly,  they  are  as  follows. 

First,  O2-  adsorbs  on  the  reactive  carbon  surface  (CRS),  and  forms 
a  strongly  bound  CO  functional  group  at  reactive  carbon  sites 
(CRS0). 

O2-  adsorption  on  the  reactive  carbon  surface  and  electron 
discharge: 

Crs  +  O2-  -  CRS0  +  2e~  (R16) 

Then,  O2-  adsorbs  at  — CrsO  sites,  and  CO2  is  produced  electro- 
chemically  (R17  and  R18).  If  the  CrsO  is  desorbed  before  the 
adsorption  of  the  second  O2-,  then  CO  is  produced  (via  R19)  [29]. 
This  group  decomposes  very  slowly  to  form  CO  (R19)  [17]. 

O2-  adsorption  at  CrsO  sites  (RDS): 

CrS0  +  O2-  -  CrsOI"  (R17) 

CO2  desorption  and  electron  discharge: 


Fig.  7.  Effect  of  bubbling  on  constant  current  discharge  curves  at  20  mA  cm  2  and 

lA/C/carbonate  =  TO  Wt/o. 


CO  desorption  (partial  electrochemical  oxidation): 


(RIB) 


CO  from  the  slurry  before  it  is  consumed  through  R12.  The  products 
of  R12  are  not  gases;  thus,  the  CO  production  rate  is  increased  by 
inhibiting  R12,  as  discussed  below. 

3.4.  Electrochemical  oxidation  of  carbon,  and  coulombic  efficiency 

Elemental  reaction  pathways  for  complete  and  partial  electro¬ 
chemical  oxidations  of  carbon  are  illustrated  in  Fig.  11,  referring 
previous  studies  [3,5,15,27,28].  The  reaction  is  initiated  by  the 
dissociation  of  carbonate  to  form  O2-  through  R14  [3,5  .  Carbon  and 
carbonate  also  react,  producing  CO  and  O2-  through  R15  [27,28  . 

CO^-^C02  +  02"  (R14) 
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Fig.  8.  Effect  of  bubbling  on  constant  current  discharge  curves  at  20  mA-2  cm  2  and 

Wc/carbonate  =  3.0  Wt/o. 


CRS0  -  CO  (g)  (R19) 

The  C/O2-  ratio  is  a  key  to  inhibit  partial  electrochemical 
oxidation.  As  noted  previously,  bubbling  provides  both  C  and  O2- 
ions  to  the  Au  surface  by  stirring,  and  it  removes  a  small  amount  of 
gas  from  the  carbon/carbonate  slurry.  When  the  C/O2-  ratio  is  low, 
as  at  WC/c arbonate  =  1.0  wt%,  complete  electrochemical  oxidation 
(R16  through  R18)  becomes  dominant,  and  Ar  bubbling  inhibits 
partial  electrochemical  oxidation  (R19),  as  shown  in  Fig.  9.  This  is 
because  O2-  ion  transport  is  dominant.  However,  when  the  C/O2- 
ratio  is  high,  as  at  Wc/C arbonate  =  3.0  wt%,  partial  electrochemical 
oxidation  (R19)  becomes  dominant,  as  shown  in  Fig.  10.  The 
bubbling  effect  becomes  insignificant,  because  carbon  transport  is 
more  dominant,  and  the  amount  of  CO  that  forms  is  beyond  the 
capacity  of  Ar  bubbling.  Therefore,  the  impact  of  bubbling  is 
dictated  by  the  C/O2-  ratio. 

Combining  the  above  elemental  steps,  global  anode  reactions 
(Rl,  R4,  and  R5)  can  be  derived.  R1  is  derived  from  R14,  R16,  R17, 
and  R18.  R4  is  derived  from  R14,  R16,  and  R19.  R5  is  derived  from 
R15,  R16,  and  R19.  These  pathways  are  described  in  Fig.  11. 

Many  researchers  have  proposed  equations  to  evaluate  the  DCFC 
performance.  In  some  literature,  carbon  efficiency  is  calculated  on 
the  basis  of  the  assumption  that  carbon  is  consumed  through  Rl 
and  R6  (the  Boudouard  reaction)  [4,13,14  .  However,  in  this  study, 
CO  is  mainly  produced  from  partial  electrochemical  oxidation. 
Therefore,  in  this  study,  it  is  appropriate  to  use  coulombic  efficiency 
(rjc),  which  is  defined  as  the  ratio  of  integrated  cell  current  to  the 
current  calculated  for  4-electron  transfer  per  atom  of  carbon  15]. 
Cooper  and  Selman  suggested  a  formula  to  calculate  coulombic 
efficiency  (Eq.  (6))  on  the  basis  of  Rl  and  R4  [15,16]: 


a  +  2 

Vc  =  2(2a  +  l) 


where  a  is  the  ratio  of  C02  and  CO  production  rate. 

Our  results  show  that  CO/C02  ratio  is  quite  higher  than  unity  at 
Wc/carbonate  =  3.0  wt%.  However,  maximum  CO/C02  ratio  is  unity 
when  Rl  and  R4  are  used.  To  explain  high  CO  production  rate,  R5 
that  produces  only  CO  during  discharge  is  necessary. 
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Fig.  9.  Effect  of  bubbling  on  measured  C02  and  CO  production  rate  during  discharge  at  20  mA  cm  2  (at  WC/ca rbonate  =  1.0  wt%). 


Thus,  in  this  research,  R1  and  R5  were  used  to  calculate 
coulombic  efficiency.  When  the  equation  of  coulombic  efficiency  is 
derived,  the  following  reactions  (R20  and  R21)  are  useful. 


and  (10)).  No  other  unknown  coefficients  can  be  added.  The  for¬ 
mula  of  coulombic  efficiency  (Eq.  (8))  can  be  rewritten  to  be 
function  of  the  ratio  of  CO2  and  CO  (a)  as  follows: 


xC  +  2XCO3  3xC02  +  4xe 
yC  +  (1  /2)yCO§-  -  (3/2)yCO  +  ye~ 


(R20) 

(R21) 


Vc  = 


(4x  +  y) 

4 


(4x  +  y) 
4(x  +  y) 


M*)l 

2{'+2©} 


a  +  2 

2(2a  +  l) 


in) 


where  x  and  y  indicate  carbon  consumption  ratios  for  Rl,  and  R5, 
respectively.  The  sum  of  x  and  y  is  given  as  follows: 

x+y=l  (7) 

Coulombic  efficiency  (rjc)  is  defined  as  n/4,  where  n  is  the 
equivalent  charge  associated  with  electrochemical  oxidation  of  one 
mole  of  carbon.  77 c,  Qco2,x,ja  and  Qco,x,y  are  derived  as  follows: 


n  (4x  +  y) 
r,c  ~  4  —  4 

(8) 

n  (3  x)/Vm 

*2C02,  product,  x,y  —  ^4x+y)p 

(9) 

(3y/2)A/m 

CiCO, product, x,y  —  (4x  +  y)F 

(10) 

In  this  approach,  the  number  of  unknown  coefficients  is  two  (x 
and  y).  The  number  of  independent  equations  is  also  two  (Eqs.  (9) 


When  Eq.  (11)  is  derived,  x  +  y  =  1  (Eq.  (7))  is  used.  Eq.  (11), 
derived  from  Rl  and  R5,  is  the  same  as  Eq.  (6),  which  is  derived 
from  Rl  and  R4.  This  is  because  R5  can  be  derived  from  Rl  and  R4 
by  eliminating  C02  mathematically.  Rl  and  R5  can  be  used  when 
the  CO/C02  ratio  is  larger  than  unity.  Qco2  and  Qco  are  also  a 
function  of  a,  as  shown  by  the  following: 

_  3  IVm 

^02, a, product  ~  2(2  +  a)F  K 


Qco, a, product 


3  aIVm 
2(2  +  a)F 


(13) 


On  the  basis  of  Eqs.  (11)  through  (13),  7)Cl  Qco2,a,Product  and 
Qco, a, product  can  be  determined  from  the  measured  CO/C02  ratio 
(=a).  This  approach  is  defined  as  “a  function”  in  this  study. 

In  addition  to  Eq.  (11),  coulombic  efficiency  can  be  determined 
from  Eq.  (8)  by  minimizing  the  error  term  defined  as  the 
followings: 
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Fig.  10.  Effect  of  bubbling  on  measured  C02  and  CO  production  rate  during  discharge  at  20  mA  cm  2  (at  WC/Carbonate  =  3.0  wt%). 
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Fig.  11.  Reaction  mechanisms  for  electrochemical  oxidation  of  carbon  in  molten  carbonate  (partially  referring  [3]). 


0cO2 ,  product,  exp  “  0cO2,x  ,y,  product 

£co2+co  = - 73 - 

*<tC02,  product,  exp 

Qco,  product,  exp  “  QcO,x  ,y, product  | 

Qco, product, exp 

In  this  approach,  carbon  consumption  ratios  (x  and  y)  that 
provide  minimum  error  term  are  determined  by  x-sweeping,  which 
is  similar  to  the  least  square  method.  Then,  Qco2,  Qco,  and  are 
determined  by  Eqs.  (8)  through  (10).  This  approach  is  defined  as 
“min.  error”  in  this  study. 

Fig.  12  shows  measured  and  calculated  CO  and  CO2  production 
rates  at  Wc/Carbonate  =  1.0  wt%.  Without  bubbling,  the  measured  CO2 
production  rate  is  lower  than  the  calculated  values  (Fig.  12(a)). 
However,  the  difference  is  not  so  significant,  this  difference  might 
be  because  of  unknown  factors  in  partial  electrochemical  oxidation 
of  activated  carbon.  Another  possible  explanation  is  that,  L^O  and 


C02  undergo  a  recombination  reaction,  given  as  R22  which  leads  to 
a  lower  CO2  production  rate. 

Li20  +  C02  Li2C03  (R22) 

Meanwhile,  bubbling  inhibits  partial  electrochemical  oxidation, 
leading  to  low  CO  formation  (Fig.  12  (b)).  Table  3  shows  estimated 
coulombic  efficiencies  and  error  terms.  Ar  bubbling  improves 
coulombic  efficiencies  by  about  10%,  regardless  of  approach  (“a 
function”  or  “min.  error”). 

Fig.  13  shows  measured  and  calculated  CO  and  CO2  production 
rates  at  Wc/carbonate  =  3.0  wt%.  Contrary  to  the  case  at  Wq 
carbonate  =  1.0  wt%,  calculated  CO  production  rate  is  50%  higher  than 
the  experimental  rate.  This  suggests  that  the  detected  CO  and  CO2 
production  rates  are  insufficient  for  discharge  at  20  mA  cm-2, 
indicating  products  other  than  CO  and  CO2  form  during  discharge. 
R12,  in  which  the  product  is  Li2C03,  may  explain  this  difference.  If 
R12  is  processed,  gas  products  are  reduced,  while  current  is  kept 
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Fig.  12.  Experimental  and  calculated  CO  and  C02  production  rate  (Wc/carbonate  = 
1.0  wt%). 
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Fig.  13.  Experimental  and  calculated  CO  and  C02  production  rate  (WJ 
carbonate  =  3.0  wt%). 


constant.  This  provides  additional  evidence  that  R12  is  processed 
during  discharge  and  at  higher  OCVs  of  1.3— 1.4  V.  It  also  indicates 
that  coulombic  efficiency  cannot  be  determined  from  R1  and  R5  at 
Wc/carbonate  =  3.0  wt%.  Therefore,  another  formula  is  necessary  to 
estimate  coulombic  efficiency  in  this  case. 

Table  4  shows  coulombic  efficiency  (rjc),  voltage  efficiency  (tjv 0i) 
and  total  electrochemical  efficiency  (r/totai)  of  the  DCFC  at  Wq 
carbonate  =  TO  wt%.  Total  electrochemical  efficiency  is  determined  by 
the  following  equations. 


nFA£ 


n 


A  £ 


AG° 


V  total 


A  H  ^theory  ^AG°^ntheoryF^ 

^  —  Vc  x  Vvol  x  ^theory ^ 


(15) 


where  carbon  (graphite)  data  is  used  as  the  theoretical  value,  nthe0ry 
is  4  and  A £  is  cell  voltage.  The  “error  min.”  approach  is  used  to 
calculate  coulombic  efficiency.  Voltage  efficiency  is  determined 
from  the  current  discharge  curve  at  the  steady  state,  as  shown  in 
Fig.  7.  Bubbling  increases  the  total  efficiency  of  the  DCFC  from  52  % 
to  64  %. 


4.  Conclusion 

This  study  aims  to  enhance  DCFC  efficiency  by  bubbling  Ar  gas 
in  a  carbon/carbonate  slurry.  The  DCFC  developed  here,  which  is  of 
a  cartridge  type,  used  a  porous  alumina  tube  to  contain  carbon/ 
carbonate  slurry  as  fuel.  The  DCFC  was  operated  at  1073  K,  and 
current  discharge  curve  and  CO  and  CO2  production  rates  were 
measured  during  discharge  at  20  mA  cm-2.  When  the  carbon 
content  in  carbonates  (Wc/carbonate)  was  1.0  wt%,  Ar  bubbling 
decreased  the  CO/C02  production  ratio  from  0.16  to  0.051,  and 
increased  the  C02  production  rate,  which  approached  the  theo¬ 
retical  value  for  complete  electrochemical  oxidation  of  carbon. 
From  the  balance  of  these  production  rates,  CO  was  mainly  formed 
from  partial  electrochemical  oxidation  rather  than  the  Boudouard 
reaction.  It  was  shown  that  bubbling  inhibited  partial  electro¬ 
chemical  oxidations.  Moreover,  bubbling  increased  the  cell  voltage 
by  approximately  0.1  V  at  steady  state.  At  Wc/carbonate  =  3.0  wt%, 
the  CO/C02  production  ratio  without  bubbling  increased  signifi¬ 
cantly  to  8.5,  indicating  that  partial  electrochemical  oxidation  was 
dominant.  The  amount  of  CO  formed  at  Wc/carbonate  =  3.0  wt%  was 
sufficient  for  CO  and/or  the  alkali  ions  to  intervene  in  the 
discharge  process,  which  showed  a  higher  OCV  (=1.3  V). 


Table  3 

Coulombic  efficiency  of  the  DCFC  at  Wc/carbonate  =  1.0  wt%. 


Without  bubbling 

With  bubbling 

a  function 

Min.  error 

a  function 

Min.  error 

Vc  [-] 

0.82 

0.86 

0.93 

0.93 

£C02+CO  [%\ 

70 

38 

3.7 

1.9 

Table  4 

Coulombic,  voltage,  and  total  efficiencies  of  DCFC  at  Wqc arbonate  =  1.0  wt%. 


With  bubbling 

Without  bubbling 

7)c  (min.  error)  [-] 

0.93 

0.86 

Vvo\  [  — ] 

0.69 

0.60 

V  total  [  ] 

0.64 

0.52 

350 
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Therefore,  voltage  at  Wc/carbonate  =  3.0  wt%  was  higher  than  that  at 
Wqc arbonate  =  1.0  Wt%.  At  Wc/carbonate  =  3.0  Wt%,  bubbling  did  not 
influence  the  current  discharge  curve,  whereas  it  did  at  Wq 
carbonate  =  1.0  wt%.  Considering  reaction  pathways,  the  C/O2-  ratio 
was  a  key  to  explain  the  electrochemical  oxidation  of  carbon. 
When  the  C/O2-  ratio  was  low  like,  as  at  Wqc arbonate  =  1.0  wt%, 
complete  electrochemical  oxidation  became  dominant,  and 
bubbling  enhanced  complete  electrochemical  oxidation  by  stir¬ 
ring.  However,  when  the  C/O2-  ratio  was  high,  as  at  Wq 
carbonate  =  3.0  wt%,  partial  electrochemical  oxidation  became 
dominant,  and  bubbling  did  not  inhibit  partial  electrochemical 
oxidation  because  carbon  transport  became  dominant.  Bubbling 
provided  both  C  and  O2-  ions  to  the  Au  surface  by  stirring. 
However,  the  impact  of  bubbling  was  dictated  by  the  C/O2-  ratio. 
Bubbling  also  promoted  removal  of  gases  from  the  carbon/car¬ 
bonate  slurry.  When  Wc/carbonate  =  1.0  wt%,  bubbling  increased 
both  coulombic  and  voltage  efficiencies,  resulting  in  an  increase  in 
total  efficiency  from  52%  to  64%. 
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Nomenclature 

a  the  ratio  of  CO2  and  CO  production  rate  [-] 

doo2  interplanar  distance  [nm] 

E0  standard  potential  [V] 

F  Faraday  constant  [C  mol-1] 

/  current  flow  [A]  (=C  s-1) 

Id  current  flow  density  [A  cnrT2] 

La  the  layer  dimension  parallel  to  the  basal  plane  [nm] 

Lc  the  layer  dimension  perpendicular  to  the  basal  plane  [nm] 

M  metal  indicating  K,  Na  or  Li 

n  equivalent  charge  associated  with  electrochemical 

oxidation  of  one  mole  of  carbon  [-] 

Qi  flow  rate  of  species  i  ( i  =  CO  or  CO2)  [ml  min”1] 

0co2,theory  theoretical  CO2  production  rate  for  complete 
electrochemical  oxidation  of  carbon  at  20  mA 
[ml  min-1] 

Qi, product  gas  production  rate  of  species  i  during  discharge 
[ml  min-1] 

Qi,x,y, product  calculated  gas  production  ratio  of  species  i  determined 
by  x  and  y  [ml  min-1  ] 

Qi, a, productcalculated  gas  production  ratio  of  species  i  determined  by 
a  [ml  min-1] 

Vm  molar  volume  of  CO2  [ml  mol-1] 

Wc/carbonate  carbon  content  in  the  carbonates  [wt%] 


Xi  mole  fraction  of  species  i  (i  =  CO  or  CO2)  [ml  min  1] 

x ,  y  carbon  consumption  ratios  for  R1  and  R5  [-] 

Greek  symbol 

cco2+co  error  term  in  Eq.  (14)  [-] 
ri  coulombic  efficiency  [-] 

A E  cell  voltage  [V] 

AG°  Gibbs  free  energy  [kj  mol-1] 

AH0  enthalpy  [kj  mol-1] 

Subscript 

c  coulombic 

exp  experiment 

vol  voltage 
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